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Abstract 

We extend the classical Gouy- Chapman model of two planar paral- 
lel interacting double-layers, which is used as a first approximation to 
describe the force between colloidal particles, by considering the finite- 
thickness of the colloids. The formation of two additional double lay- 
ers due to this finite thickness, modifies the interaction force compared 
to the Gouy-Chapman case, in which the colloids are semi-infinite ob- 
jects. In this paper we calculate this interaction force and some other 
size-dependent properties using a mean field level of description, based on 
the Poisson-Boltzmann (PB) equation. We show that in the case of finite- 
size colloids, this equation can be set in a closed form depending on the 
geometrical parameters and on their surface charge. The corresponding 
linear (Debye-Hiickel) theory and the well-known results for semi-infinite 
colloids are recovered from this formal solution after taking appropriate 
limits. We use a density functional corresponding to the PB level of de- 
scription to show how in the case when the total colloidal charge is fixed, 
it redistribute itself on their surfaces to minimize the energy of the system 
depending on the afore mentioned parameters. We study how this charge 
relaxation affects the colloidal interactions. 

1 Introduction 

Electrostatic forces play a remarkable role in the stabilization and phase behav- 
ior of colloidal suspensions. Thus, it is not surprising that understanding the 
interaction among charged colloidal particles across an intervening electrolyte 
solution has been an important goal for colloid science since the pioneering work 
by Verwey and Overbeek [[]. Poisson-Boltzmann theory (PB) allows a mean 
field calculation of these interaction at a nonlinear level. The corresponding 
linearized description, Debye-Hiickel theory (DH), is exact in the asymptotic 
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Figure 1: a) Two parallel double-layers, b) Four double-layers as a model for 
colloidal interactions. Here the finite size of the colloids is taken into account. 

limit of infinite dilutions and low-charged colloids, when the electrostatic inter- 
actions are weak compared to the thermal energy. DH theory fails to describe 
quantitatively the thermodynamic properties of realistic colloidal suspensions, 
and nonlinear PB equation must be employed, often at the expense of numerical 
computations. 

Nonetheless, the PB equation admits analytical solution in some non-trivial 
cases. An important example is that in which two parallel infinite double-layers 
interact. In this case the geometry is dictated by two planar objects (colloids) 
which are infinite in the directions parallel to their surface and semi-infinite in 
one direction perpendicular to their surface, and carrying a surface charge o. 
The region between the semi-infinite objects is filled with symmetric electrolyte 
as illustrated in figure QJi. Thus, the force between the colloids is mediated by 
the diffuse layers surrounding their surface. PB theory allows the calculation of 
this force as well as the free energy of interaction, through the determination 
of the mean local electrostatic potential in the region filled with electrolyte ^ . 
The planar geometry is relevant for colloidal spheres in the case that the Debye 
length of the electrolyte is much smaller than the radius of the spheres (such 
that the curvature effects are unimportant) . It is also relevant for understanding 
of colloidal clay platelets and their interactions. An important aspect of these 
clay platelets is their finite thickness d which can easily be of the same order as 



the Debye length and which could therefore affect the plate-plate interactions. 

The geometry of intend in this paper is depicted in figure lb, showing the fi- 
nite thickness of the platelets and hence two inner surfaces (with charge density 
Ui) and two outer surfaces (with charge density CTo). Since the plates are com- 
pletely immersed in electrolyte, double-layers form at both sides of the colloidal 
particles. The presence of charge on the outer surfaces induces the redistribu- 
tion of ions in the electrolyte, and diffuse layers of charge appear in the vicinity 
of each region in contact with the colloids. In other words, the finite-size of 
the colloidal particles allows the formation and interaction of four double-layers 
instead of two. It is clear that in this case, the force between the colloidal par- 
ticles depends not only on their relative separation h but also on their thickness 
d. 

In the traditional model of two interacting double layers ^ , the initial studies 
supposed a fixed surface charge density or a fixed potential at the surfaces. But 
there have been several studies in which the surface charge of the colloidal planes 
is not fixed, but a charge regulation mechanism can occur [21 E] , due to an ion 
adsorption-dcsorption of the colloidal surface groups. When the finite-thickness 
of the colloids is taken into account, as in our model (sec figure QJa), there is 
another kind of charge regulation that can take place. The total charge of each 
colloid can be fixed but the amounts of charge on the inner and outer surfaces 
can vary in order to minimize the thermodynamic grand potential of the system. 
We call this mechanism "charge relaxation" . This charge relaxation also affects 
the force between the colloids. In this contribution, we study quantitatively 
these finite-size effects in the case of two equally charged colloids immersed in a 
symmetrical electrolyte, in the framework of PB theory. The present nonlinear 
study already accounts for charge renormalization effects |0 El C] . 

In section|51we specify the model and introduce a suitable density functional. 
In scction|31 the nonlinear PB equation is solved. We show that for the geometry 
of figure^, the solution of the PB equation can be set in a convenient analytical 
closed form, which depends only on the geometrical parameters h and d and on 
the total surface charge density a of the colloids. The linear DH potential is 
obtained from the formal analytical solution by taking the appropriate limits. 
In section 01 we minimize the density functional to calculate the fraction of 
charge on the inner and outer surfaces, as a function of the afore mentioned 
parameters, when charge relaxation is allowed. Results are discussed in section^ 
and conclusions are given in section O 

2 The model 

We consider four infinite parallel planes located at zi = —h/2 — d, Z2 = —h/2, 
Z3 = h/2 and Zi = h/2 + d, with surface charge densities tri, tT2, <J3, respec- 
tively. The regions zi < z < Z2 correspond to the interior of the left colloid, and 
Z3 < z < Z4 to the interior of the right colloid. The regions outside the colloids 
arc filled with an electrolyte solution of dielectric constant e, at temperature T . 
Inside the colloids there is no electrolyte present but there is a medium with the 
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same dielectric constant as the electrolyte solution to avoid electrostatic image 
effects. The system is assumed to be in contact with a salt reservoir of salt den- 
sity Ps such that the electrolyte has screening length Id = = (SttAbPs)"^^^, 
with Xb = / {kBTe) the Bjerrum length, where e is the elementary charge, 
ks Boltzmann constant, and e the electric permittivity of the solvent. In what 
follows, all lengths will be conveniently given in units of the Debye length, 
namely D — Kd, H = nh and x ^ kz. For simplicity we consider the case when 
ai = (74 = (To and 02 = 03= Oi as illustrated in figure ^p. Note that with 
these conditions the system possess mirror symmetry with respect to the plane 
at .T = sketched with dashed lines. 

The thermodynamic and statistical properties of the system are character- 
ized by the dimensionless grand potential density functional w [p] = Kf3n[p]/ (2Aps) 
where A is the area of the planes and is the grand potential density functional. 
uj[p] is given by 

uj[p] = u;,d[p] + I J dxp(x)0(.T)+2/(xlogx+(l-x)log(l-x))- (1) 

The first term in the right hand side is the ideal gas entropic contribution of 
the ions in the solution, 

where p+ (x) and p- (x) are the number density profiles of the positive and nega- 
tive salt ions respectively. The second term of accounts for the electrostatic 
energy, with (j)(x) = e%l){x)/kBT the reduced electrostatic potential {jp is the 
electrostatic potential), and the total number charge in the system defined as 

4 

P{x) = P+{x) - p-{x) + -^Oj5{x ~Xj). (3) 

The last term in JQ) is of entropic origin and takes into account the fact that the 
surface charge ions on the colloids arc indistinguishable. For later convenience 
we have introduced x = ct^ /tr which is the fraction of charge in the inner surface 
with respect to the total charge density on a colloid o = Oi + Oo- Finally, 
the dimensionless constant / = AttXbo / [en) is the ratio of the Debye length 
Id = 1/k and the Gouy-Chapman length A = e/ {4ttoXb)- The later is the 
decay length of the counterions distribution around an infinite charged plane 
with bare surface charge density o in the salt free limit — > . Note that / is 
an estimate of the electrostatic energy compared to the characteristic energy of 
thermal fluctuations in the system. 

Minimization of the density functional with respect to the ion density 
profiles p±{x) gives the equilibrium density profiles in terms of the average local 
potential. In the mean field approach the reduced local electric potential (t>{x) 
is the solution of the nonlinear Poisson-Boltzmann equation 

— ^ = sinh (j) in the outside (4) 
dx'' 
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(j) 

-—r =0 in the inside (5) 
with the boundary conditions that (/> is a continuous function of x and 

0'(f')-<^'(f^)=/« (6a) 



c^\{^+D)-)-<t,'{{^^+D)+) = U (6b) 
Urn (j)'{x) = (6c) 

|a:| — ^oo 

with fi ~ A-kIbcti/ [en), fo = At:Ib<Jo/ {gk) as explained before. 

The functional |^ evaluated for the equilibrium profiles gives the equilibrium 
dimensionless grand potential. Further minimization of the later quantity with 
respect to the parameter x = Ci/cr allows to find the fraction of charge on each 
surface, when charge relaxation is allowed. This will be discussed in section 0] 
The DH level of description can be obtained by expanding to second order 
in the density and linearizing for |0| <C 1. It can also be obtained from the 
nonlinear solution as will be discussed next. 



3 Solution of Poisson-Boltzmann equation 
3.1 The Nonlinear PB Equation 

Due to the mirror symmetry it is sufficient to consider a; > since the potential 
is an even function of x: (/>(a;) = (f>{—x). Also, without loss of generality, we 
consider positively charged colloids i.e. /; > 0. The solution of the PB equation 
in the region —H/2 < x < H/2 with the boundary conditions 0(0) = (j>o and 
</i'(0) = is well-known. There are several different but equivalent expressions 
for this solution ^ (SJ |S1 |2I . After the usual first integration of the PB equation 
and using the boundary condition at x = we have 

0'(x)2^4(^sinh2^-sinh^|^) . (7) 

From this equation we see that sinh^ > sinh^ thus \(f){x)\ > |0o|- We 
deduce that (j){x) does not change its sign. Furthermore with the choice fi > 0, 
we have 4>{x) > 0oj and then <f)'{x) > 0. 

Integrating once again with respect to x and defining t ~ sinh(0/2), the 
integral can be set in the form 

.sinh(0(x)/2) 

X = / . (8) 

isinh(0o/2) Jl^t2 + l)(i2 _ sinh2(<?io/2)) 
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One recognizes an elliptic integral of the first kind JOIi defined as 



and hence we have 



1 



cosh(0o/2) 
An inversion yields 

sinh ■ 



F{9,k)= [ 



F arccos 



da 



sm a 



sinh(0o/2) 1 
sinh(0(x)/2) ' cosh((/«o/2) 



sinh(^o/2) 



cn(a;cosh(0o/2), k) ' 



(9) 



(10) 



(11) 



where cn is the Jacobian cosine amplitude. The modulus k = 1/ cosh(0o/2) of 
the elliptic functions will not be explicitly indicated in the following to lighten 
the notation. To take into account the finite-size of the colloids, we match 
this expression with the solutions in the regions H/2 < x < H/2 + D and 
H/2 + D <x. We get 



= 01 (x) = 2 sinh ^ 



sinh(0o/2) 
i{x cosli(0o/2)) 



for < a; < H/2 (12) 



{x)=Mx)=a[x^^]+MH/2) 



for H/2 <x < H/2 + D (13) 



with the dimcnsionlcss electric field inside the colloids given by 



2 sinh(0o/2) sc 



H ,00 
— cosh — 
2 2 



For the outside region x > H/2 + D wc find 

1 + Ce-=^ 



(x) =(/.3(x) = 21n 



1-Ce- 



for H/2 + D<x, 



(14) 



(15) 



with 



C; = tanh(02(f +£')/4)ef 



+ D 



(16) 

where sc(u, k) = sn(w, k) / cn(u, k) with sn(w, k) the Jacobian sine amplitude jlUL 
I12| . For a fixed 0o equations (|12|l . H13|l and (|15() give the exact analytic form 
of the mean field within the Poisson-Boltzmann theory. To completely specify 
the spatial dependence of 0, wc still need to determine its value at the origin 
00. For this wc use the boundary conditions © to get 



fo + f^ 



= sinh 



^+sinh-^ 
2 Un(fcosh^) 



sinh(0o/2) sc ( — cosh — 



2 
(17) 
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with a given by eq. . 

For given values of the parameters /i, /o, H and equation (|17() can be 
easily solved numerically for c/jq In combination with (|12|l to H16|l it gives 
the exact mean field as a function of fi, fo, and the geometrical parameters H 
and D. 

An interesting feature of 4>i{x) as given in H12|) is that if wc consider the 
prolongation of the function (j>i in cq. (|12() . for any value of x, from the general 
properties of the elliptic function cn, we notice that (f>i is a periodic function 
of X with period AK, where K = F{tt/2, k). For a; > 0, it has a first pole at 
the first positive zero of cn, at a; = if/(cosh(0o/2)). Physically, the potential 
should be bounded. This means that for a given separation H, 0o satisfies 

2K 

cosh(0o/2) 

The case H = i/max corresponds to a situation where the surface charge cr^ — > 
+00. Eq. Ipl^ serves as a guide to choose the initial value for (j)o in the numerical 
resolution of (|17|l . or in a direct numerical integration of Poisson-Boltzmann 
equation 

3.2 The Linear Regime: Debye Hiickel Theory 

We can recover the results from the linear Debye-Hiickel theory by taking the 
limit 00 ^ 1- In that limit the modulus of the elliptic functions fc — > 1. Using 
cn(M, 1) = 1/ coshu, we have (f>i{x) = 0o coshx for 0o <^ 1, which is the solution 
of the linear problem in the region < x < H/2. Similarly, using sn(7i, 1) = 
tanliM we have 

a - 00 sinh(i//2) - /„ (19) 

and equation ()17|l reduces to a linear equation which can be easily solved for 
00- The result is 

, fo+{l+D)fr 

- e"/^ + DshMH/2y 

The last equation in combination with H19() and the expression for 0i (x) , deter- 
mines also the potential in the region H/2 < x < H/2 + D, namely, (l>2{x) = 
a{x - H/2) + 4>i{H/2). Finally, hnearizing and (HHJ we get (jjsix) = [aD + 
0i(7J/2)]e^/^+^e~^. It is an easy exercise to show that the same expressions can 
be obtained with the DH approximation taken ab initio. For that, it is enough to 
take functions of the form 0i(x) = Ae^ + Be^^ , 4>2{x) = a{x — H/2) + (f)i{H/2) 
and 03(a;) = Ce^/^+^e~^. Using the boundary conditions and the continuity 
of the potential, a system of coupled linear equations for the constants A, _B, C 
and a is obtained. The solution coincides with the exact nonlinear solution in 
the limit of small potentials. 



7 




Figure 2: Exact and DH based equilibrium grand potential as a function of the 
separation distance H for two interacting finite-size planes. In this plot x — 0-5 
and / = 10. 



4 The Equilibrium Grand Potential 

As mentioned before, the total charge on the colloids is assumed to be fixed, 
but eventually different fractions of charge can move between the inner and 
outer surfaces. We took this into account in the density functional |^ by fixing 
the total charge on the colloids instead of the surface charge on the inner and 
outer surfaces separately, through the introduction of the charge fraction x. The 
functional ^ evaluated at the equilibrium density profile gives, for fixed Xi the 
grand potential of the system. Further minimization of the later quantity with 
respect to x allow us to find the fraction of charge on each surface that minimizes 
the equilibrium grand potential. The grand potential is obtained by inserting 
the expressions (|12f) . (|13|) and \ib\ into the density functional The explicit 
calculation involves well-known quadratures of the Jacobian elliptic functions, 
in particular 



du 



cn2(M, k) 



sc(u, k) dn(w, k) — i?(am(u, k), k) 



(21) 



where dn(u, k) is the Jacobian delta amplitude, E(x, k) the elliptic integral of 
second kind and am(u, k) the Jacobian amplitude ^1^]. After some calcula- 
tions we obtain 

Lu = 2(j){H/2) sinh((/)o/2) sc{H/2k) - 8 cosh((/)o/2) sc{H/2k) dn(iJ/2fc) 

+8E{a.m{H/2k),k) cosh(0o/2) - 2H smh^{(j)o/2) (22) 
+2sinh((/)(f + D)/2)^{f + D)-16 smh\<j){S- + D)/4) 
+/[2xlnx + 2(1 - x) ln(l - x) + X^iH/2) + (1 - x)4>{§ + D)]- 
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Notice that k = 1/ cosh((/)o/2), that (j){x) can be evaluated from ((T^ . ifT^ . ((TH|) . 
and that (/)o is the solution of ((T7jl and thus depends on x- For given values of 
H, D and / we determine the value of x that minimizes the grand potential H22|l 
numerically 

The DH level of description is obtained from the linear density functional 
corresponding to . In this case the dimensionless grand potential reads 

to = /[2xlnx + 2(1 - x) ln(l ~ x) + X<f>{H/2) + (1 - x)0(f + D)] (23) 

where </> is given by the expressions derived in section 3.2. Although equations 
H22(l and H23(l are a bit involved, the grand potential has the familiar repulsive 
form as that for semi-infinite colloids. This is illustrated in figure 2, where the 
linear and nonlinear results arc compared for the case / = 10 and x = 0.5. 
Some consequences of the size-dependence of the potential will be discussed in 
the next section. 

5 Results 

5.1 Charge relaxation 

Let us first consider the aforementioned effect of charge relaxation, in which 
charges can move between the colloidal surfaces. We shall refer to this as the 
annealed case, in contrast to the case in which equal amounts of charge are fixed 
("quenched") on the inner and outer surfaces. Notice that in the quenched case 
X = CTi/cr = 0.5. In figure 121 the charge fraction x is plotted as a function of 
the colloid separation H for different colloidal thicknesses. We observe that for 
inter-colloidal separations of the order of the Debye length, a fraction of charge 
moves to the outer surface. This effect substantially increases with increasing 
values of the parameter /, as we can conclude from the inset, which shows the 
fraction of charge at ^ 0. Notice that for increasing /, the values of x at 
contact tend faster to zero. Thus, the effect of increasing / is similar to what 
we observe in the main plot when D is increased, namely, a more pronounced 
relaxation of charge toward the outer surface. Notice that this effect already 
disappears at colloidal separations of the order of two times the Debye length, 
at which a small difference with the quenched case is observed. It also worth to 
mention the fact that in the main plot, no substantial change is observed when 
the thickness is increased from D = 5 to D = 10. In fact, for a given value of 
H ^ 0, when the thickness varies from D = 5 to D = oo, the charge fraction x 
does not change beyond 3% . For a colloidal thickness of a few Debye lengths, 
the distance to the outer surface is already too large to affect the fraction of 
charge that migrates as we vary the distance between the colloids. 

5.2 Electrostatic potential 

In figure^wc plot the electrostatic potential (f){x) for x > in the quenched case, 
obtained by equations H12I15|I using the methods described in section O The 
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Figure 3: Fraction of surface charge density x ^^s a function of colloid separation 
for different colloidal sizes for / = 10. In the inset we see how x{H — > 0) changes 
with varying /. 
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0.5 1 1.5 2 

X 

Figure 4: Average local electrostatic potential for H = 1, D = 0.5 and x = 0.5 
(quenched case). Different values of the dimcnsionlcss parameter / are plotted. 
Notice that the linear DH results are asymptotically exact in the limit / <C 1 
but fail to quantitatively predict the local mean field for / > 1. 
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Figure 5: Comparison of the local electrostatic potential 4>{x) for the quenched 
and annealed cases. The geometric parameters are same as in figure^] Moderate 
difference is only observed for f ^ 1 in accordance with figure 13 

dotted line represents the corresponding linear DH mean potential. Even though 
for / <C 1 both theories coincide, DH theory fails to quantitatively predict the 
local value of the mean field for / > 1. For example, for D ^ 0, (|20|l reduces to 
'Po = ifi + /o)e~^^^i thus, the total charge density ai + CTq is concentrated in a 
plane of infinitesimal width and the potential at the plane of symmetry decreases 
exponentially with respect to the plates separation. In the limit D ^ oo we 
have (j)Q = fi/ smh{H/2), this expression diverges for — » whereas eq. (|17|l 
predicts a finite value for this quantity. In figure|5l we compare the electrostatic 
field in the quenched case with that in the annealed case. Notice that as charge 
moves from the inner to the outer surface, the potential at the first decreases 
while it increases at the later. In other words relaxation of charge reduces the 
potential difference between the inner and outer surfaces. This effect is stronger 
for larger / in accordance with figure 13 

5.3 Force between the colloids 

The pressure is the force between the plates per unit area, and is defined in terms 
of the grand potential through the relation (3p/{2ps) = —du/dH. Alternatively, 
taking into account the fact that the electrostatic pressure of the system must 
balance the hydrostatic pressure, it can be written in terms of the mean field ij) 

as n 
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Figure 6: Dimensionless force per unit area between the coUoids for / = 10. 
This quantity is always positive, indicating repulsion between the colloids, and 
tends to zero for as the separation between the colloids increases. The plots for 
the quenched and annealed cases are labelled with Q and A respectively. 

where n(x) ~ P+{x) + P-{x) is the total density of microions at x. It is easy to 
show that within PB theory, this quantity gives the same value for any position 
X. In particular evaluating at x = we get 



which implies that the force between the colloids is completely determined by 
the value of the mean electrostatic potential at the plane of symmetry. This 
feature allow us to find this quantity directly from the numerical solution of 
equation (fTTI . In figureElwe plot the dimensionless pressure as a function of the 
colloidal separation, for different values of the thickness D for the quenched and 
annealed cases. Although the last two give obviously the same curve for D = 0, 
a substantial difference is observed for small, nonzero, colloidal thicknesses. On 
the other hand, for £) = 10 and larger, the inter-colloidal force turns out to be 
almost independent of the thickness for the same reasons discussed before in the 
context of charge relaxation. All the curves tend asymptotically to zero in the 
limit of infinite colloidal separation. In the DH limit (3p = 2ps(cosh0o — 1) — 
Ps(/)q. Using (PU)) wc obtain (5p ~ ps/f/sinh^ (-j) for D oo. Again, in the 
limit H 0, the linear theory predicts divergent behavior for the pressure, 
while the full nonlinear theory gives a finite value for this quantity. 

5.4 Average charge density between the colloids 

Using Poisson's equation (f)" {x) = —p{x)/{2ps), we compute the ratio (Q) of the 
average ion charge density between the two colloids, compared to the charge 



l3p = 2ps(cosh(?!)o - 1), 



(25) 
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Figure 7: Average charge density in the region between the colloids for / = 10. 
When D ^ oo charge of sign opposite to that on the colloidal surface is present 
in the middle region to neutralize the charge on the colloids. When D is finite, 
part of that charge is squeezed out of the middle region as the colloids are driven 
to a closer relative position. The insert shows the squeezing of charge when we 
decrease the thickness for fixed separation H. 



density on the inner faces of the colloids: 

(Q) = — p{x)dx = i (0'(O) - </>'(f ")) = (26) 

l^<^i Jo Ji ^ ' Ji 

This quantity is shown as a function of the inter-colloidal distance for different 
values of D in figure [7| When the distance between the finite- width colloids is 
shortened, charge can be squeezed out of the region between the colloids. This is 
not possible for semi-infinite colloids, since in that case there is no place for the 
charge to be squeezed out, such that (Q) = — 1 for all H as shown in figure [7| 
The inset shows the squeezing of charge when we decrease the thickness for 
fixed separation H. For D = the quenched and annealed cases coincide. One 
expects that Q is closer to —1 in the annealed case than in the quenched case, 
since the surface charge can follow the squeezed-out ionic charge by migration 
to the outer surface. This is indeed shown by our numerical results as shown in 
figure 7. However, for very small separations, H « 10~^ this effect is reversed, 
as can be seen in figure 8, where (Q) is closer to —1 in the quenched case. Note 
that this effect occurs for rather small values of the separation distance. In 
a model that accounts for the finite-size of the ions this effect is presumably 
absent. 
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Figure 8: Average charge density (see text) in the region between the colloids 
for / = 10. 

6 Conclusion 

The system consisting of two interacting double-layers, used as a model of col- 
loidal interactions when the curvature of the macromolecules is negligible has 
been extended to consider finite- width colloids, by the inclusion of two additional 
interacting double layers. The force that drives the colloidal particles apart is 
modified by their finite-size. The nonlinear PB equation was solved consider- 
ing this finite-size to find quantitatively the dependence of the mean field with 
respect to the geometrical parameters. We considered two relevant cases: the 
quenched case where the surface charge density on each surface of the colloids is 
the same and the annealed case, where the total charge of each colloid is fixed, 
but not the charge on each surface. In the latter case a charge "relaxation" effect 
is possible where the charges in each surface rearrange themselves in order to 
minimize the grand potential of the system. We found that this relaxation effect 
modifies also the mean electrostatic potential, and as a consequence the force 
between the colloids. We consider these results as a stepping stone towards the 
study of more complicated charge-relaxation and charge-regularization mecha- 
nisms, and to more complicated geometries, such as spheres treated within the 
Derjaguin approximation. 

7 Acknowledgments 

This work is part of the research program of the Stichting voor Fundamenteel 
Onderzoek der Materie (FOM), which is financially supported by the Neder- 
landse Organisatie voor Wetenschappelijk Onderzoek (NWO). G. T. acknowl- 



14 



edges partial financial support from COLCIENCIAS project 1204-05-13625. 

References 

[1] E. J. W. Verwey and J. Th. G. Overbeek, Theory of the stahility of lyophobic 
colloids (Elsevier, Amsterdam, 1948). 

[2] B. W. Ninham and V. A. Parscgian, J. Thcor. Biol. 31, 405 (1971). 

[3] S. H. Behrens and M. Borkovec, Phys. Rev. E 60, 7040 (1999). 

[4] S. Alexander, P. M. Chaikin, P. Grant, G. Morales, P. Pincus, and D. Hone, 
J. Chem. Phys., 80, 5776 (1984). 

[5] E. Trizac, L. Bocquet and M. Aubouy, Phys. Rev. Lett. 89, 248301 (2002). 

[6] L. Bocquet, E. Trizac, M. Aubouy, J. Chem. Phys. 117, 8138 (2002). 

[7] E. Trizac and Y. Levin, Phys. Rev. E, 69, 031403 (2004). 

[8] D. McCormack, S.L. Carnie, and D. Y. C. Chan, J. Colloid Interface 
Sci. 169, 177 (1995). 

[9] M. N. Tamashiro and H. Schiessel, Phys. Rev. E 68, 066106 (2003). 

[10] I. S. Gradshteyn, I. M. Ryzhik, Table of Integrals, Series, and Products 
(Academic Press 1994). 

[11] All the numerical resolutions (root finding and minimization) involving the 
Jacobian elliptic functions were performed using Mathematica version 5. 

[12] E. T. Whittaker and G. N. Watson, A course of modern analysis, (Cam- 
bridge, 1927) 



15 



